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مدل‌سازی فرایند پاشش حرارتی برای تشکیل پوشش سطحی با استفاده از روش‌های آماری * 
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چکیده 
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در این مقاله برای مدل‌سازی سه بعدی پوشش سطحی ناشی از فرایند پاشش حرارتی و به‌دست آوردن مشخحصات آن از رو شآماری استفاده 
شده است. در این روش, برخورد تک تک قطره‌های سرامیک مذاب و تأثی رگذار ی آن‌ها در شک لگیری لایه‌های پوشش بررسی شده است. 
برای این منظور, مدل‌های تحلیلی برای دینامیک پحش شدن قطره‌ها و برگشت لبه‌ی قطره‌های پحش شده به سمت بالابه‌کا رگرفته شد, و 
فرض شد که برگشت لبه‌های اسپلت تنها عامل ایجاد تخحلخل باشد. مدل‌سازی برای بخش کوچکی از صفحه انجام شد که به‌روش پلاسما با 
سرامیک زی رکونیای پایدارشده با یوتریوم (182) پوشش داده شده بود. در نهایت, مقدار متوسط تخلحل پوشش و ضخامت پوشش محاسبه 
شدند», و نتایج به‌دست آمده با نتایج تجربی مطایقت داشت. افزون بر اين, تأثیر عوامل محتاف فرایند شامل فاصله‌ی مشعل از سطحء سرعت 


حرکت مشعل و میزان حجم مواد مصرفی در واحد زمان بر حواص پوثشش مانند تخلحل و ضخامت نیز بررسی شده است. 


واژه‌های کلیدی مدل آماری. فرایند پافش حرارتی. تخلخل پوشش, ضخامت پوشش, برگشت لبه. 
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* نسخه‌ی نخست مقاله در تاریخ ٩۰/۳/۲۵‏ و نسخه‌ی پایانی آن در تاریخ ۹۰/۸/۸ به دفتر نشریه رسیده است. 
(۱) نویسنده‌ی مسوول: دانشجوی دکترای مکانیک. دانشگاه فردوسی مشهد 

(۲) دانشیار. گروه مهندسی مکانیک. دانشگاه فردوسی مشهد 

(۳) دانشیار. گروه مهندسی مکانیک. دانشگاه فردوسی مشهد 


۳۲ 


مقد مه 

فرایند لابه‌نشانی حرارتی یکی از روش‌های ایجاد 
پوشش سطحی با استفاده از پودر فلزات آلیاژها و 
شکل پودر و به‌وسیله‌ی روش‌های حرارتی پا انفجاری 
به ذرات مذاب يا نیمه مذاب تبدیل می‌شود (شکل ۱). 
سپس این ذرات پس از پاشش بر روی سطح مورد نظر 
به آن چسبیده و به‌این ترتیب» پوشش سطحی به‌صورت 
لایه‌هایی بر روی هم تشکیل می‌شود (شکل ۲). از این 
روش پوشش‌دهی برای تقویت و افزایش مقاومت سطح 
در برابر خوردگی و سایش. و ایجاد عایق حرارتی 
استفاده می‌شود [1]. 


شکل۱ نمایی از روش پاشش حرارتی برای تشکیل پوشش 


سطحی 


خواص مهم یک پوشش عبارتند از ضخامت» 
زبری و میزان تخلخل. این ویژگی‌ها وابستگی زیادی به 
شرایط فرایند پوشش‌دهی دارند. ذره هنگام برخورد با 
سطح می‌تواند مذاب. نیمه مذاب و یا جامد باشد 
بنابراین کیفیّت پوشش حاصل به دینامیک برخورد این 
ذره با سطح و میزان پخش شدن آن بستگی دارد. 

نتایج تجربی نشان می‌دهند که پوشش تشکیل شده 
کاملاً یک‌پارچه نبوده و بین قطره‌های پخش شده 
تخلخل وجود دارد (شکل ۲). تخلخل‌هایی که هنگام 
انجماد قطره‌ها بر روی سطح و یا لایه‌های زیرین ایجاد 


می‌شوند. به‌میزان قابل توجهی خواص مکانیکی و 


شبیه‌سازی فرآیند پاشش حرارنی جهت تشکیل پوشش.. 


تشکیل تخلخل دخالت دارند که عامل برگشتگی لبه‌ها 
به‌دلیل وجود تنش‌های حرارتی از مهم‌ترین آنها است 
[1]. این پدیده در تحقیق حاضر به‌عنوان تنها دلیل ایجاد 
تخلخل در نظر گرفته می‌شود. سازوکارهای دیگر 


سوراخ‌های موجود در لایه‌ی زیرین [2]. 


تحقیقات آزمایشگاهی زیادی در زمینه‌ی برحورد 
تک قطره‌ها به سطح انجام شده است. با در نظر گرفتن 
تأثیر عوامل مختلف مانند سرعت برخورد. دمای سطح 
و اندازه‌ی قطره بر میزان پخعش شدگی آن, رابطه‌هایی 
برای تعیین بیش‌ترین میزان پخش شدن قطره به‌دست 
آمده است [3-9]. در این رابطه‌هء ضریب پخش قطره 
( مین > ع) ارائه شده است که در آنه عم بیشن‌ترین 
قطر پخش شدگی ذره و ۲ قطر اولیه‌ی آن است. 

به‌دلیل پیچیدگی سازوکار برگشت لبه‌ی قطره‌ی 
پخش شده تحقیقات کمتری در مورد آن انجام شده 
است. فوکانوماء [10]» مدلی برای شکل گیری تخلخل 
تاش از بر کیت لبه ازایه کردهو نشان عاده: است کته 
تخلخل بیش‌تر در پیرامون قطره‌ها شکل می‌گیرد. 
به‌گونه‌ای که محل آن از حدود ۰/۶ شعاع قطره شروع 
می‌شود. سوبولف و گیلمنی [11]. یک دسته رابطه‌ی 
تحلیلی برای تشریح توزیع فشار در زیر قطره‌ی پخش 


نشریه مهندسی متالورژی و مواد 


شده به‌دست آورده‌اند که می‌توان از آن‌ها برای 
محاسبه‌ی نقطه‌ی آغاز برگشت لبه‌ها استفاده کرد. زو و 
همکارانش, [12] با مطالعه‌ی چندین مقطعم برش 
خورده از قطره‌ی پخش شده در آزمایش‌های مختلف 
میزان برگشت لبه‌ی قطره را بر حسب عوامل برخحورد 
قطره و خواص مواد به‌دست آورده‌اند. مبنای مطالعات 
آن‌ها محاسبه‌ی توزیع تنش در قطره‌ی پخش شده بوده 
است. 

ماار و همکارانش, [113 نقش عوامل مختلف موثر بر 
تخلخل موجود در پوشش تولید شده ببه‌روش 
زیرکونیای پایدار شده با یتریم را با استفاده از روش 
آماری طراحی 


آزمایش‌ها بررسی کرده‌اند. بخشی از این 


مدل‌سازی استفاده شده است. 

شناسایی عوامل موژثر بر خواص پوشش و 
چگونگی کنترل این عوامل برای به‌دست آوردن پوشش 
مطلوب. به‌کمک مدل‌سازی تشکیل پوشش با یک روش 
آماری امکان‌پذیر است. بلاش‌چنکو و چرنیاک [14] و 
سیرولینی و همکاران. [15]» از روش آماری دو بعدی 
بدون در نظر گرفتن برگشت لبه‌ی قطرات پخش شده به 
سمت بالا برای بهینه‌سازی فرایند پاشش حرارتی 
استفاده کردند. امّا نتایج آن‌ها به‌دلیل دو بعدی بودن 
مدل و استفاده از نظریه‌ی قدیمی پخش شدن قطره و 
برگشت لبه از دقّت بالایی برخوردار نبودند. غفوری آذر 
و همکاران» [16] روش آماری مونت کارلو را برای 
مدل‌سازی تشکیل پوشش بر روی سطح به‌کار بردند و 
در توزیع آماری مورد نظر برای قطره‌ها. برخورد مستقیم 
قطره‌ها به سطح را فرض کرده و زاویه‌ی برخورد قطره 
را در نظر نگرفته بودند. افزون بر این آن‌ها از توزیع 
نرمال برای زاویه‌ی برگشت لبه استفاده کرده بودند. زو 
و همکارانش. [12,17]» در ادامه‌ی کار غفوری آذر» از 
رابطه‌ی مناسب‌تری که در این تحقیق نیز از آن استفاده 
شده است. برای محاسبه‌ی زاوبه‌ی برگشت لبه‌ها 


استفاده کرده‌اند. در رابطه‌ی مذکور. خواص ترموفیزیکی 
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ماده و تأثیر اندازه‌ی قطره هنگام برخورد در نظر گرفته 
شده‌اند. امّا در مدل‌سازی از برخورد مایل قطره به سطح 
صرف نظر کرده و تنها برخورد عمودی قطره به سطح 
در نظر شده آسنتت: پنابراین» قطره‌ها بهوصورت دایره‌ای 
پخش شده و به شکل بیضی در نمی‌آیند. 

در این تحقیق. ازفک فتلال آمارش: مه بعلعن 
استفاده شده است که می‌تواند تغییرات میزان تخلخل. 
ضخامت متوسط و زبری سطح را تخمین بزند. در این 
مدل, برای هریک از متغیرهای فرایند پاشش از توزیع 
مناسب آماری استفاده شده است» و برای رفع نقعص 
موجود در کارهای قبلی افزون بر در نظر گرفتن 
برخورد مایل قطره به سطح. روش زو و همکاران 
پایان نتایج به‌دست آمده با نتایج مان [13] مقایسه 
شده است. بر اساس نتایج ارائه شده مربوط به 
حالت‌هایی که دمای قطره‌ها در پاشش کم‌تر از دمای 
ذوب است. میزان تخلخل محاسبه شده به‌دلیل وجود 
ذرات ذوب نشده کم‌تر از میزان به‌دست آمده توسط 
صحت نتایج مدل‌سازی. تأثیر عوامل مهم فرایند پاشش 
بررسی شد. ۳ ین حالت‌ها؛ یکی از عوامل 
مواد پاشیده شده و یا زاویه‌ی برگشت لبه به‌عنوان متغیر 


فرض شده و سایر عوامل ثابت در نظر گرفته می‌شوند. 


روش حل عددی فرایند پاشش حرارتی 
مدل پخش شدن قطره. در این تحقیق. همانند مدل زو 
[17) و غفوری آذر. [16]» از هم پاشیدگی قطره هنگام 
برخورد با سطح به‌دلیل رسیدن دمابه نزدیک دمای 
ذوب. در نظر گرفته نشده است. در این مدل فرض 
شده است که هر ذره قبل از برخورد با سطح به‌شکل 


کره‌ای به قطر « باشد. و بعد از برخورد به یک لایه‌ی 
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استوانه‌ای به قطر .سل تبدیل شود به‌گونه‌ای که لبه‌های 
آن به‌دلیل وجود تنش‌های حرارتی به سمت بالا برگشته 
است (شکل ۲. 

از آن‌جا که برخورد بعضی از قطره‌هابه سطح 
مایل است و پس از برخورد به‌شکل بیضی در می‌آیند. 
برای محاسبه‌ی بیش‌ترین ضریب پخش, «« هه از دو 
رابطه استفاده شده است. برای محاسبه‌ی قطر کوچک‌تر 
از رابطه‌ای استفاده شده است که توسط پسندیده فرد و 
همکاران. [18], با حل عددی برخحورد یک قطره به 
سطح و انجام آزمایش‌های مرتبط به‌دست آمده است 
(رابطه‌ی ۱- الف). برای محاسبه‌ی قطر بزرگ بیضی نیز 
از رابطه‌ی دیگری استفاده شده است که توسط اسدی و 
همکاران [20] ارائه و برای مدل‌سازی برخورد مایبل 
قطره به‌کار رفته است (رابطه‌ی ۱- ب). این رابطه‌ها با 
برقراری حالت توازن بین انرژی کل قطره قبل از 
برشورد و پس از بخش شدن آن به‌دست آمده‌اند: پسن 
از برخورد با سطح. انرژی جنبشی و انرژی سطحی 
(ناشی از تنش سطحی) اولیّه‌ی قطره برای غلبه بر اثرات 
لزجت هنگام پخش شدن آن به‌مصرف می‌رسد. 
برقراری حالت توازن بین انرژی جنبشی. کشش سطحی 
و انرژی تلف شده نرخ پخش شدن قطره را تعیین 
می‌کند. این رابطه به‌وسیله‌ی محققّان زیادی استفاده شده 


شبیه‌سازی فرآیند پاشش حرارنی جهت تشکیل پوشش.. 


همکاران [16) ستقیمی و همکاران» [21 و عابدیشی 
و همکاران» [22) اشاره کرد. این رابطه‌ها به‌شکل زیر 
۳ 


(۱-الف) برخورد عمودی قطره 
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در این رابطه‌ها؛ اعداد بدون بعد ری‌نولدز (6) 
وبر (۷۷۵). پرانتل (:۳). استفان )٩16(‏ و پکلت (۳6۵) به 


وسیله‌ی رابطه‌های زیر تعریف می‌شوند: 
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قطره‌ی پخش شده 
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قطره قبل از برخورد 
زاویه‌ی برگشت لبه 
سطح چسپیده 
به دیواره‌ی قطره 


شکل۳ تصویر ذره قبل و بعد از برخورد. و نمایش برگشت لبه‌ی قطره (: محل شروع برگشت لبه. و ۵ : زاویه‌ی برگشت لبه) 
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افزون بر این 1 زاویه‌ی برخورد قطره با سطم, 0 
زاویه‌ی تماس بین مایم و سطح تحت پوشش؛ 
شک( رب سای ولهق سظم: عضیریب: مایت 
حرارتی. ۷ ضریب کشش سطحی, ۸ لزجت. 6 گرمای 
ویژه‌ی مذاب 1 دمای ذوب و :13 گرمای نهان ذوب 


ذره هلر 


مدل شکل گیری لایه 

در این تحقیق. چگونگی قرارگیری لایه‌های ناشی از 
پخش ذرات بر روی هم به‌کمک کامپیوتر و بر مبنای 
نت خضمین شبیال ۴( ۱۲۵۲ سل سازین شاه اس 
[16,17]. در این روش از تابع اسکالر ۶ برای هر سلول 
متاسیاین اسفاقه نله شک مه کور مدای که کنر ستلول 
پر از سبال باشده مقذان #برابر با یک خواهد برد مقدار 
برای سلول‌های خالی برابر با صفر و بر روی مرز 
قطره در بازه‌ی ۱ >۶ > 0 قرار دارد (شکل ۴). 


مدل تحلیلی برگشت لبه‌ی اسپلت 


لبه‌ی قطره‌های پخش شده در نتیجه‌ی وجود تنش‌های 


پس‌ماند حرارتی به‌سمت بالا بر می‌گردد. به‌دلیل وجود 
اختلاف دما بین دیواره و قطره. قطره بعد از پیخش 


۳۵ 


شدن سرد شده و شروع به انجماد می‌کند. با توجه به 
تفاوت نرخ انجماد در لایه‌ی نزدیک به دیواره و لایه‌ی 
نزدیک به سطح آزاد و با در نظر گرفتن وجود ضریب 
انقباض حرارتی ماده در نتیجه‌ی انجماد. تنش‌های 
پس‌ماند در اسپلت ظاهر می‌شوند. این پدیده یکی از 
مهم ترین دلایل به‌وجود آمدن تخلخل در پوشش است. 
دو عامل اساسی یعنی محل شروع و زاویه‌ی برگشت 
آن» در برگشت لبه‌ی اسپلت دخالت دارند (شکل ۳. 
محل شروع بالا رفتگی و میزان آن, به جنس ذره و 
شرایط برخورد. و نیز به عواملی مانند اختلاف ضریب 
انبساط حرارتی فاز مایع و جامد. کشش سطحی فاز 
مایع» زبری سطح و ذوب شدن مجدد ذره بستگی 
دارند [3]. مشاهدات فوکانوما. [10] نشان می‌دهد که 
اکثر حفره‌ها در کناره‌ی اسپلت و در فاصله‌ی حدود ۶۰ 
همکاران» [16]. از همین مقدار در تحلیل خود استفاده 
کرده‌اند. ولی آن‌ها برای محاسبه‌ی زاوبه‌ی برگشت 
لبه‌ی قطره از رابطه‌ی خاصی استفاده نکرده و به جای 


آن, از توزیع آماری ویژه‌ای استفاده کرده‌اند. 


(الف) 


(ب) 


شکل ۴ الف) برگشت لبه‌ی قطره در محاسبات به‌روش نسبت حجمی سیال (۷۵۵۲). ب) روش نسبت حجمی سیال برای پیدا کردن سطح 


ازاد قطره 


م4 


شکال این روش حل. فرض نادرست علدم وابستگی 
ژاو تفع بر کیت یه فتعامت: اسات انیت هدر انسماه 
آن تأثیر زیادی دارد. افزون بر این» غفوری آذر نتایج 
خود را با نتایج آزمایشگاهی مقایسه نکرده و تنها بر 


ستفاده از روش‌های آماری تأکید داشته است. با این 


حال کار او سرآغاز خوبی برای مدل‌سازی بوده است. 
زو و همکاران» [12]» برای بهبود کار غفوری آذر 
رابطه‌ی ساده‌ی زیر را برای تعیین زاویه‌ی برگشت لبه 
به‌دست آورده‌اند که در این تحقیق نیز از آن استفاده 
شده است: 


(۳) كت ك < 0 


5 


در این رابطه. ,9 ضخامت اسپلت. 0 ضریب 
انبساط حرارتی, 8/2 شعاع اسپلت. ۸ اختلاف دما 
بین قطره و سطح و « نقطه‌ی شروع برگشت است 
(شکل(۳)» که همگی بر اساس داده‌ههای تجربی 
به‌دست آمده‌اند. در تحقیق زو و همکاران, [12] صئه 
گذاری مدل با مقایسه‌ی نتایج آن با نتایج تجربی 
حاصل از بررسی مستقیم تخلخل, انجام شده است. 


مدل آماری 


سه فرض اساسی در یک مدل آماری عبارتند از: عدم 
تداخل ذرات مذاب با یکدیگر تصادفی بودن پاشش 
ذرات (هر ذره دارای اندازه سرعت ۷ و محل برخورد 
مشخحص است). و احتمال وجود یک ذره در هر لحظه 
مستقل از سایر ذرات. با اين‌که پاشش فرایندی وابسته 
به زمان است و طبیعتاً باید از تعریف زنجیره‌ی مارگف 
برای تعیین متغیرهای تصادفی مربوط به پارامترهای 
فرایند استفاده کرد امٌا فرایند برای ساده‌سازی 
یکنواخت فرض شده است. با توجه به شرایط فوق, در 
این مدل فرض شده است که سرعت. دمای قطره و 
زاویه‌ی پاشش از تابع توزیع احتمال نرمال (رابطه‌ی۴) 
و اندازه‌ی قطره‌ها از تابع توزیع احتمال لگاریتمی 
نرمال (رابطه‌ی ۵) پیروی کنند [۱۷ و ۱۶]. 


شبیه‌سازی ف رآیند پاشش حرارتی جهت تشکیل پوشش.. 


(۴) تابع توزیع نرمال 


کج اره له 
2 27 

(۵) تابع توزیع لگاریتمی نرمال 

یت ره ِ < (),ع 
20 2( 
در اين رابطه‌ها؛ ع و بو متفیرهای احتمالی» # 
میانگین» » انحراف معیار ‏ لگاریتم انحراف معیار و 
میانگین لگاریتمی هستند. توزیع نرمال به‌شکل یک 
نافوّش,می‌باشف که حول میزان میانگین آن (:2) متقارن 
است. میزان پهنی این ناقوس با انحراف معیار نشان 
داده می‌شود. بسیاری از پدیده های طبیعی به‌طور 
تقریبی از توزیع نرمال پیروی می‌کنند. توزیع لگاریتمی 
نرمال در آمار و احتمال توزیعی است که لگاریتم 
طبیعی آن دارای توزیع نرمال با پارامترهای 2 وه 
می‌باشد. به‌عبارت دیگر اگر 26 متغیری با توزیع نرمال 
باشد. (65006 < ۲ دارای توزیع لگاریتمی نرمال 
خواهد بود. 

پودر مورد استفاده در مدل‌سازی از جنس 
زیرکونیای پایدار شده با یوتريم بود. این پودر با 
علامت تجاری 204015 1۷6/60 ٩126۲‏ و مشخصات 
ص16 ع< مره ص50 ع< وق و ص95 < موله دارای 
ذراتی به‌شکل کروی بود. با توجه به اين که مشخصه‌ی 
95 < موه به معنی آن است که ٩۰‏ درصد ذرات 
پودر قطری کم‌تر از 0 ۹۵ دارند. و هم‌چنین سایر 
مشخصات پودر توزیع آماری مورد نظر را می‌توان 
تعیین کرد. افزون بر این ذرات پودر با عنوان 
1۲ 125+1- 27028۷260 مشخص می‌شوند. و معنی 
آن این است که درون پودر ذراتی به قطر کوچک‌تر از 
سب ۱۱ و بزرگ‌تر از ۰« ۱۲۵ وجود ندارند. بنابراین 
در مدل‌سازی این تحقیق پس از به‌دست آوردن قطر 
ذرات. ذراتی را که قطر آن‌ها بین ۱۱ و ۱۲۵ میکرون 
هستند با استفاده از تابع آماری لگاریتمی نرمال (مشابه 


شکل 6۵ انتخاب می‌کنيم. 
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روش حل 
پس از شبکه‌بندی میدان حل, ابتدا با استفاده از روش 
آماری عوامل مورد نیاز مشخص می‌شوند. برای پیدا 
کردن قطر اسپلت از رابطه‌ی ۲ استفاده کرده و با 
به‌کارگیری روش ۷0۳ مرز اسپلت در میدان حل 
پس از پخش شدن قطره با توجه به 


رابطه‌ی ۳ و فرض کردن مقدار ثابت ۰,۶ برای 4 


تعیین می‌شود. رد 


زاوبه‌ی برگشت محاسبه شده و از آن در جابجایی 
مرزهای اسپلت استفاده می‌شود. این عمل منجر به 
خالی شدن فضای زیر لبه‌های آن می‌شود. این کار برای 
قطره‌ی بعدی هم انجام می‌شود. با این تفاوت که 
هنگام استفاده از روش ۷0 برای پیدا کردن مرز 
قطره. وجود قطره‌های قبلی درون میدان حل در نظر 
گرفته می‌شود و پس از پیدا شدن موقغیت سطح 
زیرین» چگونگی پخش شدن قطره‌ی جدید بررسی 
هن موه 

برای تولید اعداد تصادفی. از الگوریتم ارائه شده 
ترتنط الن میا اسفاده کنتم استه گن تهاست پس از 
پایان مرحله‌ی پاشش قطره‌هاو برای محاسبه‌ی 
ضخامت و میزان تخلخل. الگوریتم‌های متفاوتی به‌کار 
برده شدند. رابطه‌ی زیر برای تعیین میزان تخلخل 


استفاده شد؛ 


۷ 
(ع 00 شش - نز005 
۷ + ,۷ ِ ۳ 


ِ 10 


که در آن ۷ فضای خالی محاسبه شده بین قطره‌های 
پضش ستسالن ۵ و ۷ حجم اشغال تخل ۵ توسط ماده 


ی بان 


نتایج و بحث 
شکل (۶) آنچه که هنگام پوشش دهی اتفاق 
می‌افتد را در چند مرحله نشان می‌دهد. در این شکل. 
پراکندگی آماری برخورد قطره‌ها با سطح مشاهده 


می‌شود. همان گونه که مشاهده می‌شود. برخورد ۵۰۰ 


۳۷ 


قطره‌ی آلومینا با سرعت متوسط لته 100 و قطر 
متوسط ۵۰ میکرون در مراحل مختلف. پس از برخورد 
۰ ۲۰۰ ۳۰۰و ۵۰۰ قطره با سطح نشان داده شده 
است. 
صحخه گذاری تسایج. برای صحه گذاری نتایج 
مدل‌سازی, از داده‌ها و نتایج ماار و همکاران» [۱۳ 
استفاده شد. در این شرایط مشخصات فرایند اسپری 
پلاسما (۸۵۳5) برای زیرکونیای پایدار شده بایوتريم 
(92) گو نطر گرقشه شسته ایسن واه از جسسن 
5 ۱۲۵/۱6 9۷1767 (با مشخصات ۸۸ 16 < 410 
۶ 57 < 50 و ۸۰ 95 < 090) دارای توزیع لگاریتمی 
نرمال با میانگین حل ۵۷ و انحراف معیار ۱/۴۷ بوده‌اند. 
شرایط انجام آزمایش به‌همراه نتایج آن در جداول (۱) 
و(۲) آمده‌اند. 

همان‌طور که در نتایج حدول (۲) مشاهده 
می‌شود. مدل‌سازی انجام شده در این تحقیق از 
توان‌مندی خوبی برخوردار بوده است. در این مدل. اگر 
چه میزان تخلل محاسبه شده کم‌تر از مقادیر محاسبه 
شده توسط ماار می‌باشد اشا روند تغیبرات تفریبا 
بکسان: استت: 

افزون بر اين» برای اطمینان از توان‌مندی مدل 
به‌کار رفته نتایج زو و همکاران نیز بررسی شدند و 


نتایج این بررسی در حدول (۳) نشان داده شده است. 


۱۱5۱09۲۵۳ 0۲ 1 
۱۱۵۲۵۱ 
3000 1 ۳ 
0 
1 132170 


مت وتا 
۳ 
دح 
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ی هم ۲ 
ی کی ی ی یی نیقی 
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شکل۵ توزیع لگاریتمی نرمال مربوط به اندازه‌ی ذرات 


۵ 393 34 عووع»اهز) جح ۴ 75 46.7 ۵992 ۱6۷۴() ۱۷2 


۴ 44۱ 54و وه مماه زو »1/2 ان ۳ 
نت 95 8 هل 
شکل ۶ طرح‌واره‌ای از پوشش دهی در سطحی به ابعاد 1۳۳0150 
(قطر متوسط ذرات ۱0 50 و سرعت برخورد متوسط ذرات با سطح 10/5 1000 بوده‌اند) 
جدول ۱ مشخصات فرایند پوشش سطح 
سرعت پاشش پودر صنص/ع۲۵ 
فاصله مشعل از سطح متغیر از ۲۰۰ تا ۲۵۰ میلیمتر 
سرعت حرکت مشعل 1/۹ 2( 
هسیک سر کیت سقیها | دنت وی کت 
جدول ۲ مقایسه‌ی با نتایج ماار و همکاران [13] 
فاصله‌ی مشعل از سرعت تخلخل (از ضخامت (۸۲2) ضخامت 
۱ دمای ذرات ۱ تخلخل 
ردیف سطح ۵ ذرات نتایج ماار و دای (از نتایج ماار و (حصب) 
(رصصص) (۳/5) همکاران) همکاران) (محاسبه شده) 
۱ ۳.۰ ۳۷۶۰ ۱۸۹ ۹/۳۳ ۹/۲ ۳۸ ۴۳۰۰ 
۳ ۳.۰ ۱۳۸۴ ۱۹۷ ۸/۳ ۸/۱ ۴۳۰۷ ۳۸۷ 
۳ ۳۵۰ ۳۷۰۲ ۱۸۳ ۹/۳۳ ۹۸۳ ۳۸ ۳۰۳ 
جدول ۳ مقایسه‌ی نتایج با نتایج زو و همکاران [17] 
۱ دمای ذرات سرعت ذرات قطر ذرات تفت ات ِِ 
تِ 0 له ۳ ۳ ۳( ۰ 
زو و همکاران) و همکاران) با رابطه‌ی ۶) 
۱ ۳۷۶ ۱۶۸ ۶۱ ۶ ۱/۸۷ ۱/۵ 
۳۲ ۳۱/۳۹۴ ۱۸۱ 0۰ ۸/۸ ۷ و2 
۳ ۱۳۸۱ ۳۰۶ ۳۹ 2/۸ 8/۶ 0/۱ 
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همان گونه که مشاهده می‌شود در اين مقایسه نیز 
میزان تخلخل محاسبه شده از رابطه‌ی ۶ کم‌تر از میزان 
آن در نتایج زو است. امّا دارای روند یک‌سانی بوده و 
در مقایسه با نتایج محاسبات زو به واقعیت نزدیک 
است. 
بررسی عوامل مهسم فرایند پاشش. پس از کسب 
اطمینان از صحّت مدل‌سازی صورت گرفته. به بررسی 
تأثیر عوامل مختلف فرایند بر میزان تخلخل. ضخامت 
و زبری سطح پوشش داده شده پرداخته‌ایم و نتایج 
حاصل با نتایج به‌دست آمده توسط زو مقایسه شده‌اند. 
فاصله‌ی مشعل از سطح عاملی است که فرایند مورد 
نظر را تحت تاثیر خود قرار می‌دهد (شکل‌های (۷) و 


100 150 200 


۳۹ 


(۸)). در مدل‌سازی‌ها. سرعت مشعل بتزآنو نیب 
5 ۰ و دبی جرمی پودر برابر با 21/5 ۰/۱۳ 
فرض شده‌اند. شکل (۷) نشان می‌دهد که فاصله‌ی 
مشعل از سطح تاثیر زیادی بر میزان تخلخل نداشته و 
در آن تخلخل تقریباً ابت مانده است. در این حالت 
نیز تخلخل محاسبه شده کم‌تر از میزان آن در نتایج زو 
بوده است. 

در شکل (۸) دیده می‌شود که با افزايش فاصله‌ی 
مشعل. ضخامت کاهش می‌یابد. دلیل کاهش ضخامت. 
برخورد کم‌تر ذرات به محدوده‌ی مورد نظر از سطح 


می‌باشد. در این حالت. مطابقت روند تغییرات با نتایج 


زو چشم‌گیر است. 


۴۵۲۵5۲۰۷ )۶۵( 
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شکل ۸ تغییرات تخلخل بر حسب فاصله‌ی مشعل از سطح 


(سرعت مشعل برابر با 10505 400 دبی جرمی پودر برابر با 81/5 0.13 و قطر متوسط ذرات برابر با تنل 25) 
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برای نمایش تاثیر محاسبه‌ی زاویه‌ی برگشت لبه 
بر محاسبات در شکل‌های )٩(‏ و (۱۰ زاویه‌ی 
کشت لبه اسظنانا از زانط (۳ ماس نله سیف 
اگر عه این خالت واقغی نسته اما برای تشان داذن 
تأثیر زاویه و نحوه‌ی محاسبه‌ی آن ضروری است. به 
این ترتیب. برای زاویه‌های مختلف با انحراف معیار 
ثابت ۵ و فاصله‌ی ثابت مشعل از سطح برابر با 2۳ 
۰ میزان‌های تخلخل و ضخامت محاسبه شده‌اند. 


4 5 6 7 8 


در شکل )٩(‏ نشان داده شده است که زاویه‌ی پرگشت 
تأثیر زیادی بر ضخامت محاسبه شده ندارد اما از شکل 
(۱۰) می‌توان دریافت که زاویه‌ی برگشت لبه تأثیر 
بسیار زیادی بر تخلخل محاسبه شده دارد. با افزایش 
زاویه‌ی برگشت لبه. میزان تخلخل به‌دلیل ایجاد فضای 
خالی زیاد در زیر اسپلت. افزایش زیادی می‌یابد. 
بنابراین؛ انتخاب رابطه‌ی مناسبی برای محاسبه‌ی این 


زاویه از اهمئّت ویژه‌ای برخوردار امتتتا: 


() ۲۳016۱0656 
زر 
‌ 


0 1 2 


(۵۲66ع06) عاع۸ ولا عطنا۲نت 


شکا ٩‏ تغییرات ضخامت بر حسب زاویه‌ی بررگشت لبه 
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(۱68۲۶6) اه ولا عصناست 


ذرات برابر با تلم 25 و فاصله‌ی مشعل از سطح برابر با هه 150) 


نشریه مهندسی متالورژی و مواد 
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(فاصله‌ی مشعل از سطح برابر با 10107 150 دبی جرمی پودر برابر با 8/5 4.13 قطر متوسط ذرات برابر با 1 25 و زاویه‌ی متوسط 
برگشت لبه برابر با 4) 
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لبه پرابر با *۵) 


برای بررسی تأثیر سرعت حرکت مشعل. آنرا بین ۱/۳ 
تا ۰/۶ متر بر ثانیه تغییر داده و برای جلوگیری از تأثیر 
تغییرات زاویه‌ی برگشت لبه در این مرحله. آن‌را ثاببت 
و برابر با ۴ درجه فرض کردیم. نتایج به‌دست آمده با 
حفظ فاصله‌ی ۱۵۰ میلی‌متر مشعل از سطم در 
شکل‌های (۱۱) و (۱۲) ارائه شده‌اند. 

همان‌گونه که در شکل (۱۱) دیده می‌شود با 
افزایش سرعت مشعل ضخامت کاهش می‌يابد. دلیل 
منطقی آن را می‌توان برخورد تعداد کم‌تری قطره در 


زمان ثابت پاشش دانست که منجر به کم‌تر شدن میزان 
ماده‌ای می‌شود که بر روی سطح می‌نشیند. از طرف 
دیگر. در شکل (۱۲) دیده می‌شود که سرعت حرکت 
مشعل تأثیری بر میزان تخلخل محاسبه شده ندارد. 
دلیل عدم تأثیر میزان تخلخل» سرعت نسبی بسیار کم 
حرکت مشعل نسبت به سرعت برخورد قطره‌هابا 
سطح است که باعث ثابت ماندن زاویه‌ی برخورد قطره 
با سطح می‌شود. به‌این ترتیب. روند مطابق انتظار است. 

و شک ای( ۱۲ و( ۱۳ فان تن دی 


۳۲ 


جرمی پودر پاشیده شده به سطح بر عوامل فیزیکی 
پوشتن مان خاقم ده اسیه: یل ( ۱۲ ان اس دهل 
که با افزایش تبین جرهی پنودر» سامت افزایش 
می‌یابد. دلیل آن را می‌توان. همانند شکل (۱۱ مقدار 
جرم موادی دانست که بر روی سطح می‌نشیند. واضح 
اس گهبا آقرایش فین جرف وان .مراد باشیکه شتده 
بر روی سطح بیش‌تر می‌شود و در نتبجه. ضخامت 
افزایش می‌یابد. 

شکل (۱۴) نشان می‌دهد که با تغییر دبی. میزان 
تخلخل ثابت می‌ماند. از آن‌جا که دبی تأثیری بر 
برگشت لبه ندارده بر روی فضای خالی زیر لبه‌های 


اسپلت و در نتیجه, بر میزان تخلخل نیز تأثیری ندارد. 


تسیر ری 20 


15 "۰-75 


شبیه‌سازی فرآیند پاشش حرارتی جهت تشکیل پوشش.. 


نتیجه گیری 
در این پژوهش. فرایند پاشش حرارتی به‌روش آماری 
مدل‌سازی شد و در آن, تأثیر عوامل مختلف فرایند بر 
روی تخلخل ب‌دست آمده و ضخامت پوشش بررسی 
شد و نتایج زیر ن4دسق: آملنكه 
- مقایسه‌ی نتایج حاصل از این تحقیق که با 
داده‌های تجربی موجود و نتایج عددی تحقیق 
دیگر در جدول (۲) صحه گذاری شده است 
ی خی که کی ختوون فتتی تفس 


برای مدل‌سازی پوشش‌دهی سطح می‌باشد. 
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شکل ۱۳ تغییرات ضخامت بر حسب دبی جرمی پودر 


(فاصله‌ی مشعا از ابر با 09۳0 150 سرعت مد اي با 90/5 400 قطر متو سط ذرات براير با فطل 25 زاویه‌ی متو سط بر گشت 
ی ر‌ اقا ی انیت او 2 مب ر‌ زاین ۶ ۳ زاویه‌ی متو بر 


لبه برابر با *4) 
20 
5 _ 
کوج 
ب 
سسسسسست 0 2 
۵ 
5 
6۵ 
5 
0 
00175 0015 0025 0013 00۰75 
(8۳/5) ۲۵۲ ۶660 ۳۴۵۱۷۷16۲ 


شکل ۱۴ تغییرات ضخامت بر حسب دبی جرمی پودر 


(فاصله مشعل از سطح-15010۳0: سرعت مشعل< 7۳/5 400 قطر متوسط ذرات> 25۳0 زاویه متوسط برگشتگی لبه<4۴) 


نشریه مهندسی متالورژی و مواد ّ 


۲- نتایج به‌دست آمده از این تحقیق نشان می‌دهند تاثیری بر میزان تخلخل ندارند. 
که نحوه‌ی محاسبه‌ی زاویه‌ی برگشت لبه تأثیر ۴ افزون بر این نتایج نشان دادند که با افزایش 
زیادی بر نتایج دارد (شکل (۱۰)). سرعت مشعل و فاصله‌ی مشعل از سطح. 
۳- نتایج نشان می‌دهند که عواملی از قبیل فاصله‌ی ضخامت پوشش کاهش می‌پابد. 


مشعل از سطح. سرعت مشعل و دبی جرمی پودر 
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